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COASTAL  ACOUSTIC  TOMOGRAPHY  DATA 
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Abstract 

A  direct  insertion  scheme  for  assimilating  coastal  acoustic  tomo- 
graphic (CAT)  vertical  temperature  sections  into  a  multilevel,  coastal 
primitive  equation  model  for  the  Gulf  of  Sirte  (Lybia)  is  investigated 
using  computer  simulation  experiments.  Although  the  model  was 
developed  for  the  whole  Mediterranean  Sea  (MED),  only  a  model 
sub-domain  covering  the  Gulf  of  Sirte  was  used  in  this  study.  The 
model  has  realistic  coastlines  and  bottom  topography,  and  a  coastal- 
following,  curvilinear,  nearly  orthogonal,  horizontal  coordinate  system 
with  a  horizontal  resolution  of  about  10  by  10  km.  The  grid  of  the 
model  was  designed  using  a  grid  generation /focusing  technique.  The 
model  has  complete  thermodynamics,  second  order  turbulence  closure, 
and  16  bottom-following  (sigma)  vertical  levels.  To  generate  the  "true 
ocean"  for  this  study,  the  model  was  first  spun  up  for  30  days  with 
the  Levitus  temperature  and  salinity  and  ECMWF  wind  climatolo- 
gies, and  then  run  for  one  year  and  more.  The  last  60  days  of  this 
control  run  were  taken  to  represent  "actuality." 

A  series  of  assimilation  experiments  was  carried  out  in  which  CAT 
temperature  sUces  synthesized  from  different  CAT  configurations  based 
on  the  "true  ocean"  were  inserted  into  the  model  at  various  time  steps 
to  examine  the  convergence  of  this  direct  insertion  scheme.  In  all  the 
assimilation  experiments,  after  a  spinup  for  30  days,  the  model  was 
then  integrated  with  CAT  slices  inserted  daily  to  produce  "nowcast 
fields"  for  periods  of  30  and  60  days. 

Our  results  strongly  indicate  that  the  coastal  model  and  the  direct 
insertion  scheme  for  CAT  slices  work  well  together.    As  time  pro- 


gresses,  the  CAT  slices  spread  information  out  in  nearly  all  directions 
in  the  Gulf  of  Sirte.  The  inserted  data  act  as  an  initial  convergence 
maker  which  forces  the  model  to  "lose"  memory  of  the  initial  state 
and  gradually  converge  to  the  "true  ocean".  The  global  rms  error  of 
the  nowcast  temperature  fields  in  the  sub-  domain  considered  in  an 
experiment  involving  five  CAT  sections  decreases  by  about  50%  after 
30  days  and  66%  after  60  days  of  data  assimilation.  The  rms  errors 
of  two  other  experiments  involving  only  a  single  slice  are  reduced  by 
33%  after  30  days.  In  all  three  experiments,  the  rate  of  convergence, 
as  measured  by  the  rms  error,  is  approximately  constant.  It  is  ex- 
pected that  the  convergence  to  the  "true  ocean"  will  be  quicker  with 
more  frequent  insertions  and  with  more  CAT  sections. 

It  is  seen  from  these  experiments  that  the  temperature  information 
spreads  out  faster  over  the  shallow-water  areas.  This  may  be  due 
to  stronger  mixing,  horizontal  shear  and  interaction  with  topography 
resulting  in  an  increase  of  horizontal  diflPusion.  The  temperatures  at 
the  mixed-layer  level  converge  more  rapidly  to  the  "true  temperatures" 
in  comparison  to  the  surface-level  temperatures.  All  the  CAT  slices 
have  almost  the  same  effectiveness  in  improving  the  nowcast  fields 
regardless  of  their  horizontal  orientations.  Our  experiments  also  show 
that  convergence  is  not  linearly  proportional  to  the  number  of  CAT 
slices  assimilated.  Our  simulation  results  confirm  that  the  internal 
forcing  provided  by  vertical  slices  derived  by  CAT  is  not  localized  in 
physical  space.  They  induce  immediate  correlation  between  the  fully 
nonlinear  primitive  equation  model  and  the  "true  ocean". 


1     INTRODUCTION 

In  the  past  ocean  modeling  tended  to  proceed  separately  from  observations, 
mainly  due  to  the  lack  of  synoptic  observations  to  properly  constrain  the 
model  physics.  Two  new  technologies  in  today's  Oceanography  brings  new 
large  synoptic  datasets.  Satellites  provide  global  maps  of  sea  surface  height 
(SSH),  sea  surface  temperature  (SST)  and  wind  stress.  The  latter  field  is  the 
most  important  surface  forcing  of  the  ocean  circulation.  The  satellites  give 
information  on  the  physical  characteristics  of  the  ocean  surface  and  the  near- 
surface  layer.  On  the  other  hand,  Ocean  acoustic  tomography  can  provide 
quality  4-D  data  measuring  the  interior  ocean  structure.  This  underwater 
acoustic  inverse  technique  uses  travel  time  changes  of  sound  pulses  to  map 
sound  speed/temperature  perturbation  along  various  paths.  The  satellites 
and  ocean  tomography  are  highly  complimentary  ocean  observing  systems. 
The  ability  to  efficiently  assimilate  these  datasets  into  ocean  models  is  key 
to  ocean  nowcasting  and  forecasting. 

Ocean  tomography  has  traditionally  been  used  in  deep  water  for  mapping 
the  "ocean  weather."  The  adaptation  of  this  technique  from  deep  to  shallow 
water  poses  many  scientific  and  engineering  challenges.  New  sensor  systems, 
advanced  sound  propagation  models,  and  optimal  signal  processing  and  map- 
ping techniques  that  will  work  in  complex  coastal  environments  must  be  de- 
veloped and  tested  at  sea.  Research  into  and  development  of  coastal  tomog- 
raphy is  the  main  thrust  of  a  collaboration  between  scientists  and  engineers 
at  the  Naval  Postgraduate  School  (NPS)  and  Woods  Hole  Oceanographic 
Institution  (WHOI).  In  August  1992,  a  first-generation  coastal  acoustic  to- 


mography  (CAT)  system  was  tested  in  the  Barents  Sea  (Chiu  et  al.,  1993, 
and  Miller  et  al.,  1993).  This  system  consists  of  an  electronically  controlled 
sound  source  and  a  telemetered  vertical  hydrophone  array  (Von  Der  Heydt 
et  al.,  1992).  Using  such  system,  the  two  groups  demonstrated  that  CAT  is 
practical  and  can  produce  vertical  slices  of  the  coastal  ocean  temperature  and 
sound  speed  fields  at  high  temporal  and  spatial  resolutions.  An  assessment  of 
the  potential  performance  of  a  coastal  nowcasting  system  which  assimilates 
CAT  vertical  slices  into  a  state-of-the-art  coastal  ocean  model  is  the  focal 
point  of  a  computer  simulation  study  reported  here. 

In  the  assimilation  of  oceanographic  datasets,  numerical  ocean  models 
play  a  critical  role  in  the  advection  and  propagation  of  the  information  con- 
tained in  the  data.  Two  outstanding  questions  were  raised  by  Ghil  and 
Malanotte-Rizzoli  (1991):  (1)  Can  information  provided  only  at  the  sea  sur- 
face be  transferred  dynamically  into  the  deep  oceanic  layers,  thus  reconstruct- 
ing the  deep  circulation?  (2)  Can  information  provided  only  locally,  in  limited 
oceanic  regions,  be  transferred  to  ocean  areas  far  away  from  the  data-dense 
region  and  which  time  and  space  scales  are  better  estimated  through  the  as- 
similation of  local  data?  These  issues,  particularly  pertaining  to  deep-ocean 
circulations,  have  been  addressed  by  Ghil  and  Malanotte-Rizzoli  (1991)  as 
well  as  other  experts  working  in  the  subject  area.  We  will  give  a  brief  review 
of  some  of  these  previous  computer  simulation  works  which  are  somewhat 
related  to  our  present  work  in  studying  the  direct  insertion  approach  and 
advection  information  process.  The  central  difference  between  this  work  and 
the  previous  works  is  that  this  work  focuses  on  nowcasting  in  the  coastal 
oceans  using  CAT  and  a  state-of-the-art  coastal  ocean  model. 


The  problem  of  propagation  of  surface  information  to  the  ocean  inte- 
rior was  first  studied  by  Hurlburt  (1986)  and  Thompson  (1986).  They  both 
inserted  the  altimeter  data  directly  at  every  grid  points  into  a  two-layer  prim- 
itive equation  models  for  the  Gulf  of  Mexico.  Focusing  on  the  dynamic  trans- 
fer of  information  from  the  surface  to  the  deep  layer,  Hurlburt  (1986)  demon- 
strated the  success  of  the  numerical  ocean  model  used  in  the  reconstruction 
of  deep  circulation.  Hurlburt  showed  that  approximately  two  updates  per 
eddy  cycle  (  57  days)  are  required  for  the  assimilation  to  be  successful  and 
provide  convergence  to  the  reference  ocean  and  the  convergence  is  clearly 
accelerated  further  when  updating  every  20  days.  Thompson  (1986)  focused 
on  the  geoid  error  as  it  affects  the  assimilation  of  altimeter  data  for  mesoscale 
ocean  prediction.  He  showed  that  errors  in  the  geoid  on  spatial  scales  com- 
parable to  the  model  gird  resolution  did  not  seriously  degrade  the  forecast, 
even  in  dynamically  active  regions  with  large  gradients  in  the  geoid  height 
where  instability  process  has  been  observed  to  occur.  Kindle  (1986)  used  a 
one-layer  reduced-gravity  primitive  equation  model  of  the  Gulf  of  Mexico  to 
examine  the  methodology  of  incorporating  satellite  altimeter  data  into  ocean 
forecasting  model.  His  main  result  was  that  for  a  stationary  circular  eddy, 
approximately  two  track  (either  ascending  or  descending)  across  the  eddy 
are  sufficient  to  ensure  adequate  spatial  resolution  and  an  irregularly  shaped 
eddy  may  require  three  or  four  tracks. 

In  a  two-part  study,  Malanotte-Rizzoli  and  Holland  (1986  and  1988;  here- 
after MRH)  examined  the  transfer  of  information  by  advection  from  data  lo- 
calized in  space  to  other  regions  of  the  ocean  using  an  idealized,  mid-latitude, 
multilayer,  quasi-geostrophic  (QG),  eddy- resolving  model  and  a  direct  inser- 


tion  technique,  weighting  observations  by  their  distance  from  the  grid  point 
being  updated.  In  the  first  part  (MRH,  1986)  they  used  a  quasi-linear,  steady 
state  model  ocean  and  the  assimilated  data  were  density  measured  along  vari- 
ous hydrographic  sections.  The  main  result  of  the  study  is  that  a  local  section 
can  be  quite  effective  in  determining  flow  in  regions  far  away  from  the  data 
section  if  the  flow  is  simple,  steady,  and  quasi-linear.  The  most  effective  data 
sections  are  meridional,  long  and  far  away  from  the  ocean's  western  bound- 
ary. In  the  second  part  (MRH,  1988),  they  used  a  time-dependent,  fully 
eddy-  resolving  model  ocean  and  the  assimilated  data  were  sections  with  dif- 
ferent orientations.  They  found  that,  on  the  short  time  scale  of  mesoscale 
variability,  all  the  data  sections  considered  are  equally  ineffective.  A  single 
data  section  is  quite  effective  in  driving  the  model  to  the  reference  ocean  if 
the  data  insertion  process  is  carried  out  for  time  longer  than  the  model  equi- 
libration time.  A  single  section  of  assimilated  data  is  completely  ineffective 
in  driving  the  model  toward  reality  over  time  scales  comparable  with  the 
mesoscale  variability.  Their  result  also  show  that  the  model  climatological 
mean  after  ten  years  of  data  assimilation  becomes  extremely  similar  to  the 
climatological  mean  of  the  reference  ocean. 

For  this  present  study,  we  use  a  highly-nonlinear,  primitive  equation, 
eddy-resolving,  coastal  model  for  the  Gulf  of  Sirte  (the  Gulf  of  Lybia). 
The  model  uses  curvilinear  coordinates  with  coastal-following  and  nearly- 
orthogonal  grid.  The  model  physics  include  fully  active  thermodynamics 
and  a  second  order  turbulence  closure.  This  coastal  model  was  developed  by 
adapting  the  Princeton  Ocean  Model  (POM)  for  the  Gulf  Stream  (Mellor  and 
Ezer,  1991)  to  the  Mediterranean  Sea  (MED).  Previously,  Mellor  and  Ezer 


(1991)  used  POM  and  an  optimal-interpolation  altimeter-data  assimilation 
scheme  to  study  the  quality  of  the  nowcast  and  forecast  fields  in  the  Gulf 
Stream  region.  Here,  we  use  the  model  to  study  the  potential  of  using  CAT 
for  coastal  ocean  nowcasting.  Our  assessment  is  focused  on  the  simple  direct 
insertion  approach.  An  interesting  question  is  that  what  configurations  of 
CAT  in  a  shallow  water  region  such  as  the  Gulf  of  Sirte  (Gulf  of  Lybia), 
where  both  topographic  and  nonlinear  effects  are  strong,  are  most  effective, 
i.e.,  lead  to  a  large  and  widespread  improvement  of  the  nowcast  fields.  Our 
emphasis  here  is  on  the  understanding  of  horizontal  propagation  of  the  infor- 
mation obtained  by  CAT,  the  effect  of  data  insertion  intervals  on  the  nowcast 
fields,  and  the  convergence  rate  associated  with  the  direct  insertion  scheme 
using  this  state-of-the-art  curvilinear,  multilevel  primitive  equation  coastal 
ocean  model  for  the  Gulf  of  Sirte. 

The  report  is  organized  as  follows.  In  the  next  section,  we  give  a  descrip- 
tion of  coastal  acoustic  tomography  (CAT)  with  a  focus  on  the  inversion  of 
tomographic  data  to  obtain  vertical  temperature  slices.  The  coastal  ocean 
circulation  model  for  the  Gulf  of  Sirte  is  described  in  Section  3.  The  model 
was  first  used  to  simulate  the  "true  ocean"  and,  subsequently,  it  was  used  to 
produce  nowcast  temperature  fields  with  the  assimilation  of  synthetic  CAT 
slices  for  various  CAT  configurations.  The  direct  insertion  technique  is  also 
discussed  in  Section  3.  In  Section  4,  the  numerical  experiments  and  the 
results  are  discussed.  Section  5  gives  a  summary  of  our  findings. 


2     COASTAL  ACOUSTIC  TOMOGRAPHY 

In  contrast  to  deep  water  propagation,  the  arrivals  from  the  multipaths  due 
to  a  pulse  excitation  in  a  shallow-water  environment  tend  to  overlap  in  time. 
It  is  therefore  difficult  to  resolve  them  individually  with  omni-directional 
or  almost  omni-directional  receivers  which  are  typically  used  in  deep  water 
tomographic  applications.  Vertical  line  arrays  and  beamforming  techniques 
can  be  used  to  overcome  this  difficulty.  The  observed  perturbations  of  the 
beamformed  individual  ray/mode  travel  times  constitute  the  database  of  the 
coastal  tomographic  inverse  scheme  developed  by  Chiu  et  al.  (1993).  Their 
scheme  is  discussed  next  in  detail. 

2.1      A  Hybrid  Ray/Mode  Inverse  Method 

Establishing  the  mathematical  relations  between  the  observations  and  the 
unknown  structure  is  the  first  step  in  the  development  of  any  inverse  scheme. 
In  our  case,  the  observations  are  the  changes  in  ray  and  mode  travel  times, 
6Vj  and  6i^,  and  the  unknown  ocean  variable  is  perturbation  of  sound  speed, 
6c,  i.e.,  the  deviation  from  a  reference  sound  speed  field  c.  The  reference 
ocean  is  allowed  to  be  range  dependent  in  our  formulation. 

2.1.1     Ray  Travel  Time 

The  relation  of  the  travel  time  of  a  pulse  signal  along  a  raypath,  i^,  to  the 
speed  of  sound  c  is  well  known.  The  relation  is 


r.=  [  -^ds  (1) 


where  x  =  (a:,  y,  z)  is  the  position  vector  and  s  is  the  arclength  along  the  path 
trajectory  Tj  connecting  a  source  and  a  receiver.  Replacing  c  by  c-\- 6c  and  t^j 
by  i^j  -f  SVj  and  assuming  the  change  in  travel  time  due  to  path  perturbation 
is  of  2nd  order,  we  obtain,  through  linearization, 

2.1.2     Mode  Travel  Time 

The  derivation  of  the  relation  between  bf^  and  8c  is  less  straight  forward.  The 
uniformly  valid,  asymptotic  coupled-mode  solution  of  Desaubies  et  al.  (1986) 
shall  be  our  starting  point.  This  solution  was  obtained  by  expanding  both 
the  magnitude  and  phase  of  the  modal  horizontal  structure  in  a  perturbation 
series  in  powers  of  a  small  parameter  //,  where  //  is  the  ratio  of  the  wavelength 
to  the  characteristic  horizontal  scale  of  the  oceanic  variability.  The  horizontal 
phase  of  the  n-the  mode,  to  0(/i'^),  is  given  as 

where  r  is  the  range  along  the  horizontal  mode  path  connecting  a  source  and 
a  receiver,  «„  and  k^  are  the  horizontal  wavenumbers  of  modes  n  and  m, 
respectively,  and  Amn  is  a  coupling  coefficient  that  can  be  equated  as 

Amn=         —. T[Zn{z;r,uj) Z^(2;r,tu) \dz    (4) 

Jo     p{r,  z)  or  or 

In  the  depth  integral  of  (4),  p  is  density,  and  Z^  and  Zm  are  the  nth  and 
mth  modes,  respectively. 

Analogous  to  the  "static  ray  approximation,"  we  assume  that  the  varia- 
tions in  the  mode  functions  produce  higher-order  changes  in  the  horizontal 
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phases  comparing  to  those  caused  directly  by  the  variations  in  the  horizontal 
wavenumbers  due  to  a  perturbation  in  sound  speed.  This  "static  mode  ap- 
proximation" was  introduced  by  Rajan  et  al.  (1987)  in  their  inverse  method 
for  obtaining  geo  parameters.  Under  this  assumption,  the  variation  of  the 
modal  phase  integral  (3)  leads  to 


80 


TOO 

n=    I       {SKu  +  Yl  fmn^^n  +  Yl  9mn^l<'m)dr  (5) 


with 


-     -       7        2 
ymn  T-    2         -      2^2  ^    ■' 

where  the  modal  quantities  with  an  "overbar"  are  those  associated  with  the 
reference  ocean.  From  a  stationary  phase  consideration,  modal  travel  time 
and  phase  are  related  by 

C  =  ^  (8) 

It  follows  from  (8)  and  (5)  that  the  variation  in  modal  travel  time  is 

dSOn 


''^    =      du 


06 Kn     ,    v-^  dfmn  .  ,    V"^   r        ^<^«n 


To  find  the  relation  between  bt"^  and  <5c,  we  need  the  relations  of  6/c„ 
and  hence     a   "  to  be.  Using  the  governing  vertical-structure  equation  and 
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orthonormal  condition  for  the  normal  modes,  Raj  an  et  al.  (1987)  showed 
that  the  perturbation  of  «„  due  to  a  small  8c  is 

8Kn  =    /        QriScdz  (10) 

with 

qn  =  — 3^-^71  (11) 

Note  that  the  "static  mode  approximation"  was  also  imposed  in  obtaining 
the  eigenvalue  perturbation  equation.  Finally,  a  substitution  of  (10)  in  (9) 
gives 

<^C     =     /     /     { ^ +  2^— [/„,„(r,2,a;)g„(r,z,u;) 

-\-gmn{r,z,uj)qm{r,z,uj)\}8c{x)dzdr  (12) 

Previously,  Shang  (1989)  has  developed  a  tomographic  inverse  technique 
based  on  adiabatic  mode  theory.  His  formula  for  modal  travel  time  per- 
turbation has  only  the  first  term  in  the  integrand.  Here,  we  have  included 
the  effects  of  mode  coupling  in  the  travel  time  change. 

2.1.3      Optimal  Inverse  Solution 

We  represent  8c  in  the  vertical  slice  traversed  by  the  rays  and  modes  as  a 
Fourier  series,  i.e., 

7/2  A72 

^c=     ^        ]^     a,j7;  exp[>/^(zXoa;  +  A:Vo2)]  (13) 

.=-7/2  k--Kl2 


with 


C  =  f  (14) 

K  =  I  (15) 
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where  X  and  Z  are  the  horizontal  and  vertical  extents  of  the  slice,  respec- 
tively, and  (7  -j-  1)  x  [K  +  1)  is  the  total  number  of  spectral  components 
used.  This  truncated  Fourier  decomposition  of  8c  is  equivalent  to  gridding 
(i.e.,  discretizing)  the  physical  space  (a:,  z)  into  (/+  1)  x  [K  +  1)  points  with 
a  horizontal  grid  spacing  of  ^  and  a  vertical  grid  spacing  of  ^.  Using  the 
Fourier  decomposition,  the  inverse  problem  can  be  cast  into  simple  linear 
algebra. 

Putting  the  Fourier  coefficients  a.^'s  in  a  column  vector,  a,  and  collecting 
all  the  measured  ray  travel  time  changes  along  resolved  multipaths  connect- 
ing sources  and  receivers  in  a  column  vector,  8f  ^  the  corresponding  data- 
unknown  relation  becomes 

8f  =  G''a-^e'  (16) 

where  G^  is  the  ray  data  kernel  matrix,  calculated  by  numerically  integrating 
(2)  for  each  of  the  Fourier  components.  Similarly,  the  data-unknown  relation 
corresponding  to  the  modes  can  be  cast  as 

^'^  =  G™a -I- e"^  (17) 

where  the  data  vector  6f^  contains  the  measured  mode  travel  time  changes 
and  the  mode  data  kernel  matrix  G"^  is  computed  from  (12).  In  general,  the 
data  are  contaminated  by  experimental  noise,  a  combination  of  measurement 
and  model  errors.  We  represent  the  contamination  in  (16)  and  (17)  by  the 
noise  vectors  e'"  and  e"*,  respectively. 

There  are  several  ways  to  derive  the  linear,  minimum  mean- square-error 
estimate,  i.e.,  inverse  solution  (Cornuelle,  1985,  and  Chiu,  1985).  Assuming 
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that  a,  e^  and  e^  have  zero  means  and  are  uncorrelated,  it  can  be  shown  that 
the  minimum  mean-square-error  solution  can  be  obtained  by  minimizing  the 
following  objective  (or  cost)  function: 

+ia'^Ca-^a  (18) 

where  Cg,  C™,  and  Ca  are  the  covariance  matrices  of  e*",  e*"  and  a,  respec- 
tively. Furthermore,  the  error  covariance  of  the  solution  is  simply  the  inverse 
of  the  Hessian  (the  matrix  of  second  derivatives)  of  S  (Chiu,  1985).  One 
can  see  that  the  objective  function  is  constructed  by  a  combination  of  two 
data  constraints  (associated  with  rays  and  modes,  respectively)  and  a  statis- 
tical information  constraint.  It  is  worth  mentioning  that  if  e*",  e*"  and  a  had 
normal  statistics,  the  minimum  of  (18)  would  coincide  with  the  mode  of  the 
a  posteriori  distribution,  which  is  the  probability  distribution  function  of  a, 
given  6f  and  St"^.  In  such  case,  the  minimum  mean-square-error  estimate 
and  the  mode  of  the  a  posteriori  distribution  estimate  would  be  identical 
(Chiu,  1985). 

The  minimum  mean-square-error  inverse  solution,  a,  and  its  error  covari- 
ance, C^,  follow  from  the  evaluations  of  the  first  and  the  second  derivatives 
of  S  with  respective  to  a.  They  are 


a    =    C,{G^'Cl-'6f-\-G'^'Cf-'Sr-)  (19) 

a    =    (C^Cr'G^-hG^^C^-'G^^-hCa-')-'  (20) 

It  is  straight  forward  to  go  back  to  the  physical  space  {x,z)  from  the 
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Fourier  domain.  The  inverse  solution  for  sound  speed  perturbation,  ^c,  and 
its  error  covariance  are  obtainable  through  applications  of  multi-dimensional 
discrete  Inverse  Fast  Fourier  Transforms  (IFFT)  to  their  Fourier-space  coun- 
terparts. Ocean  temperature  perturbation  is  linearly  proportional  to  sound 
speed  perturbation.  Thus,  the  calculation  of  temperature  estimates  from 
sound  speed  estimates  is  straight  forward. 

3     COASTAL  OCEAN  CIRCULATION  MODEL  FOR  THE  GULF 
OF  SIRTE 

3.1     The  Model 

A  coastal  primitive  equation  model  with  a  fully  active  thermodynamics  is  de- 
veloped for  whole  Mediterranean  Sea  (MED)  based  on  the  Princeton  Ocean 
Model  (POM)  for  the  Gulf  Stream  (Mellor  and  Ezer,  1991).  This  is  a 
three-dimensional  MED  model  with  a  horizontal  coastal-following,  curvilin- 
ear nearly  orthogonal  and  bottom-following,  sigma,  vertical  coordinate  sys- 
tem. The  horizontal  coastal-following,  curvilinear  nearly  orthogonal  model 
grid  is  designed  using  the  EAGLEView  software  package  developed  at  Missis- 
sippi State  University  (MSU).  While  the  package  allows  for  a  high-resolution 
handling  of  the  open  boundary,  coastlines  and  Gulf  regions,  a  reduction  of 
computational  expense  is  still  achieved  (Ly  et  aL,  1993).  The  model  has 
free  surface  dynamics  and  it  uses  a  second  order  turbulence  closure  scheme 
for  the  mixed  layer  and  a  split  mode  time  step.  The  model  has  one  open 
boundary  in  the  Atlantic  ocean  about  200  km  west  of  the  Strait  of  Gibral- 
tar.   Only  the  model  sub-domain  covering  the  Gulf  of  Sirte  is  used  in  this 

15 


data  assimilation  study.   The  average  horizontal  resolutions  for  the  studied 
region  is  approximately  10  km. 

The  model  equations  governing  the  velocity,  surface  elevation,  salinity, 
and  temperature  fields  in  the  ocean  have  been  described  by  Blumberg  and 
Mellor  (1983;  1987).  The  equations  are  written  in  a  system  of  Cartesian 
coordinates  with  x  eastward,  y  northward,  and  z  upward.  The  free  surface  is 
located  a,t  z  =  r]{x,y,t)  and  the  bottom  is  at  2:  =  —H{x^y).  The  equations 
are 


f +  (7.Vi/  +  W'f +  2nxf7  =  -lvP  +  ^(A'„f)  +  F        (21) 

at  OZ  po  OZ  OZ 


dP 

a7  =  -"*  (22) 


-      dW 
V  •  ^  +  V  =  0  (23) 

OZ 


^  +  U-Ve,  +  Vvf  =  ^iA  +  F,  (24) 

Ot  OZ         OZ  OZ 

The  density  is  computed  using  the  equation  of  state  (Mellor,  1991)  in  the 
general  form 

p  =  pe{S,9,p)  (25) 

The  Coriolis  force  is  denoted  as  2fl  x  U,  where  J7  is  the  earth  rotation 
vector,  U  is  the  horizontal  velocity  vector  with  components  {U,V),  V  is  the 


16 


horizontal  gradient  operator,  po  is  the  reference  density,  p  is  the  in  situ  den- 
sity, g  is  the  gravitational  acceleration,  P  is  the  pressure,  Kyn  and  Kh  are 
the  vertical  turbulent  exchange  coefficients  for  momentum,  and  for  heat  and 
salt,  respectively.  In  equation  (24),  Oi  may  represent  mean  potential  temper- 
ature, 9  (or  in  situ  temperature  for  shallow  water  application)  or  salinity,  S. 
The  potential  density  used  here  is  an  approximation,  since  it  excludes  the 
effects  of  pressure  variations  (Blumberg  and  Mellor,  1981).  Here,  F{Fj;,Fy) 
in  equation  (21)  and  Fg,  in  equation  (24)  are  the  horizontal  mixing  terms 
(see  Blumberg  and  Mellor,  1987). 

The  momentum  conservation  and  diffusion  equations  (21)  and  (24)  con- 
tain the  vertical  turbulent  exchange  coefficients  which  are  determined  by  a 
second  order  turbulence  closure  scheme  of  Mellor  and  Yamada  (1982).  The 
turbulence  scheme  is  characterized  by  equations  for  turbulent  kinetic  energy 
(TKE),  9^/2,  and  for  the  turbulent  mixing  length,  £.  The  two  equations  can 
be  written  in  the  same  form  (see  Ly,  1992a),  so  that  Qi  is  either  q^/2  for 
TKE  or  q^£  for  the  turbulent  mixing  length  in  the  equation: 

--^U.VQ.^W—  =  _(/^,_)+/^^^[(_)  +(_)  ]^K,^---K,^-^F^ 

(26) 
where 

KmQ  =  2Km,    Khq=2Kh,    KQ=2/e,    Fq  =  F,  (27) 

when  Qi  represents  [q^  12) 
and 

KmQ=iE^Km,    KhQ=eEJ<,,    Kq  =  W,    FQ  =  Fe  (28) 
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when  Qi  represents  {q^^).  Also  in  the  above  equations,  VV^  is  a  "wall  proxim- 
ity" function  (  Mellor  and  Yamada,  1982)  and  Bi  is  an  empirical  constant. 

The  first  term  on  the  right-hand  side  of  (26)  is  the  diffusion  contribution, 
the  next  two  terms  represent  shear  production,  the  next  is  the  buoyancy  term, 
the  next  term  represents  dissipation,  and  the  last  represents  the  horizontal 
mixing.  The  vertical  turbulent  exchange  coefficients  (TEC)  Km-,  I^hi  and  Kq 
are  defined  as  Km  =  ^QSm,  Kh  =  ^qSh,  Kg  =  ^qSq  where  Sm,Sh,  and  Sg 
are  the  stability  functions  (see  Mellor,  1973;  Mellor  and  Yamada,  1982). 

At  the  free  surface,  z  =  T]{x,y),  the  surface  wind  stress,  heat,  and  salinity 
fluxes  are  prescribed.  At  the  bottom  z  =  —H{x,y),  zero  heat  and  salinity 
fluxes  are  used.  Turbulent  mixing  length,  £,  vanishes  at  the  bottom.  At 
land  boundaries  the  condition  of  no  diffusive  fluxes  of  any  property  across 
the  interface  is  used.  Sigma  coordinates  can  adequately  model  domains  like 
the  Gulf  of  Sirte  with  large  bathymetric  irregularities.  Mode  splitting  of  the 
barotropic  and  baroclinic  modes  are  imposed  in  the  model  to  save  computer 
resources. 

3.2      Model  Parameters  and  Initialization 

The  Gulf  of  Sirte  off  Lybia  is  located  in  the  southern  middle  part  of  the 
Mediterranean  (MED).  The  region  of  interest  in  this  study  is  bounded  by 
latitudes  30  and  35  N  and  longitudes  15  and  22  E.  The  eastern  part  of  the 
Gulf  of  Sirte  is  a  deep  region  where  the  100  —  m  isobath  is  less  than  20 
km  from  the  coastline.  The  southern  and  western  parts  of  the  Gulf  are 
considerably  shallower.  The  bottom  depths  along  the  CAT  sections  used  for 
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the  assimilation  are  typically  from  200  m  to  2000  m. 

In  the  model,  the  Gulf  of  Sirte  and  its  surroundings  have  a  horizontal  res- 
olution of  about  10  X  10  km.  There  are  16  bottom-following  (sigma)  vertical 
levels  with  increased  resolution  in  the  mixed  layer  and  lower  resolution  in  the 
deeper  layers.  The  external  (barotropic)  mode  time  step  is  15  5,  and  the  in- 
ternal (baroclinic)  one  is  10  minutes,  satisfying  the  Courant-Friedrichs-Levy 
(CFL)  computational  stability  criterion.  The  model  bathymetry  is  obtained 
from  the  global  DBDB5  bathymetry  dataset  with  5'  x  5'  resolution  (National 
Geophysical  Data  Center,  1985)  and  is  interpolated  to  the  model  grid.  The 
bathymetry  is  filtered  with  a  Shapiro  (1970)  filter  to  remove  high  frequency 
noise.  The  model  is  initialized  with  the  Levitus  temperature  and  salinity 
climatology  (Levitus,  1982),  and  wind  stress  from  the  ECMWF  (Ly  et  a/., 
1992b)  climatology.  In  all  our  computer  simulation  experiments,  the  model  is 
spun  up  for  30  days  diagnostically  (temperatures  and  salinities  held  constant) 
before  making  the  prognostic  run. 

3.3     Data  Direct  Insertion  Technique 

There  are  various  data  assimilation  approaches.  The  simplest  and  the  most 
computationally  tractable  is  the  direct  insertion  approach  which  is  used  in 
our  investigation.  In  the  following  discussion  of  the  direct  insertion  technique, 
we  use  superscripts  /,  o,  and  a  to  denote  the  forecast,  observed,  and  analyzed 
(the  best  first  guess)  ocean  fields,  respectively.  We  also  assume  that  CAT 
data  is  available  at  time  intervals  of  At. 

Generally,  the  ocean  field  predicted  at  time  <,  prior  to  the  assimilation  of 
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the  data  obtained  at  ^,  is  related  to  the  analyzed  ocean  field  at  time  t  —  At 
by 

Tf  =  r(Tt^,)  (29) 

where  F  is  the  mathematical  operator  describing  the  ocean  dynamics,  which 
in  our  cast  is  nonlinear.  In  the  direct  insertion  method,  the  inserted  data  can 
be  viewed  as  a  "perturbation"  introduced  at  a  specific  location  {xs,ys,Zs}. 
We  insert  data  at  all  time  stepes.  The  vector  T  represents  the  discretized 
ocean  variables  such  as  velocity  and  temperature  at  the  ocean  model  grid 
points.  An  observational  vector  can  be  thought  of  as  consisting  of  a  predicted 
and  a  perturbation  part  e  as: 

T?  =  T[  +  £  (30) 

Then  in  the  assimilation  experiments  we  have 

r(Tt  At)  =  T«  =  Tf  +  £  (31) 

at  specific  locations  {x,y,z}  =  {0:3,1/3,23}.  And 

r(TtAt)  =  T[  (32) 

at  all  other  locations.  Where  £  is  a  difference  between  the  predicted  and 
observed  value  (a  perterbation)  and  grows  in  time  and  space.  It  propagates 
and  diffuses  away  from  the  CAT  slice  section.  It  is  noted  that  the  "true 
ocean"  can  be  produced  by  replacing  the  predicted  value  T[  with  an  observed 
T°  at  spatial  and  temporal  resolutions  which  are  dense  and  frequent  enough 
to  force  the  difference  e  toward  zero. 
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The  analyzed  field  at  time  t  —  At  in  the  insertion  experiments  can  be 
constructed  by  combining  the  forecast  field  and  the  observed  field  at  t  —  At, 
that  is 

Tt  At  =  (I  -  W)T{_^t  +  WT?_^t  (33) 

where  I  is  an  identity  matrix  and  W  is  a  diagonal  weighting  matrix.  The 
diagonal  elements  of  W  consist  of  zeros  and  ones.  The  element  is  one  if  the 
corresponding  model  grid  point  is  traversed  by  the  CAT  sections  and  if  the 
corresponding  variable  is  temperature.  The  element  is  zero  otherwise.  The 
construction  of  the  analyzed  field  using  equation  (33)  and  the  prediction  of 
the  forecast  field  using  equation  (29)  is  repeated  at  every  data  assimilation 
time  step. 

4     NUMERICAL  EXPERIMENTS  AND  DISCUSSIONS 

In  addition  to  the  model  run  that  generated  the  "true  ocean"  (which  will 
be  referred  as  the  TR  experiment),  four  other  model  runs  were  performed. 
The  latter  runs  will  be  referred  as  the  EXn  experiments  with  n  =  0, 1,2,3. 
EXo  is  an  experiment  with  no  data  inserted.  This  experiment  provides  the 
baseline  for  judging  improvements  when  data  are  inserted.  In  EXi  (fig- 
ure 10),  the  model  was  run  with  a  single  CAT  section  inserted  daily.  This 
CAT  section  has  a  northeast-southwest  (NE-SW)  orientation.  In  EX2  (figure 
17),  a  single  CAT  section  was  used  again  but,  this  time,  the  section  has  a 
northwest-southeast  (NW-SE)  orientation.  In  EX3,  the  model  was  run  with 
the  insertion  of  a  total  of  five  CAT  sections.  The  five  sections  include  the 
previous  two  diagonal  slices  as  well  as  three  other  slices  enclosing  a  coastal 
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volume  making  a  CAT  parallelogram  (figure  24).    These  experiments  are 
discussed  in  detail  next. 

4.1     The  "True  Ocean"  Experiment  TR 

The  scheme  of  the  twin  experiments  in  this  study  is  shown  in  Table  1.  The 
central  issue  is  the  convergence  problem  in  using  a  fully  nonlinear  primitive 
equation  coastal  model  with  the  direct  data  insertion  technique.  The  global 
volume  rms  errors  between  the  temperature  fields  of  various  experiments  and 
the  "true  ocean"  can  be  measures  of  the  convergence  in  the  domain  studied. 
The  model  is  spun  up  for  30  days  and  then  run  for  a  year  with  ECMWF 
winds  to  produce  the  "true  ocean"  TR.  One  year  is  run  to  make  sure  the 
model  "loses"  the  initial  state  (the  loss  of  memory  of  the  initial  state).  Fig- 
ures 1,  2,  and  3  show  the  "true  ocean"  temperature  fields  for  day  one  at  the 
surface  (level  1),  the  mixed  layer  level  (level  6),  and  the  bottom  level  (level 
15),  respectively.  Figures  4-6  show  the  "true  ocean"  temperature  fields  at  the 
above  levels  for  day  30.  Hereafter,  in  the  TR  run,  days  1  and  30  correspond 
to  days  361  and  390,  respectively.  After  one  year  of  simulation  the  model  is 
stable  with  the  annual  ECMWF  winds.  From  figures  1  and  4  it  is  seen  the 
SST  fields  are  very  slightly  developing  after  30  days  (after  one  year  of  model 
simulation).  The  SST  field  in  the  shallow  waters  (the  southern,  southwest- 
ern, and  southeastern  regions  of  the  Gulf)  stay  very  much  the  same  after  30 
days  of  simulation.  Temperatures  in  these  regions  are  about  13  deg  C .  In  the 
central  and  northern  parts  of  the  region,  the  SST  field  is  slightly  different. 
The  region  of  12  deg  C  has  developed  towards  the  northwest  after  30  days 
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simulation.  The  colder  water  of  the  deeper  layer  has  moved  up  to  the  surface 
by  the  mixing  process  in  the  central  and  northern  parts  (the  deepest  region 
of  the  MED  with  depths  of  about  4000  m). 

Figures  2  and  5  show  the  temperature  field  of  the  "true  ocean"  at  level  6 
for  one  and  30  days,  respectively.  It  is  noted  that  level  6  of  sigma  coordinates 
corresponds  to  depths  from  the  subsurface  to  100  m.  The  temperature  fields 
in  all  parts  of  the  region  stay  very  much  the  same  after  30  days  simulation. 
They  have  changed  little  compared  to  the  surface  level  (figs  1  and  4).  The 
temperature  fields  are  almost  identical  at  the  bottom  level  (level  15)  after  30 
days  simulation  which  are  shown  in  Figures  3  and  6. 

As  expected,  the  30-day  run  changes  slightly  the  surface  temperature 
field.  In  the  deeper  layers  (levels  6  and  15),  the  temperature  fields  stay 
almost  the  same  after  the  30-day  run. 

4.2      The  No  Data  Assimilated  Experiment  EXo 

The  model  is  spun  up  for  30  days  with  Levitus  temperature  and  salinity  and 
ECMWF  wind  climatologies.  The  no  data  assimilated  experiment  EXg  is 
performed  by  running  the  model  for  30  days  without  the  CAT  data  insertion. 
The  temperature  fields  without  data  insertion  at  days  one  and  30  are  shown 
in  Figures  7  and  8,  respectively.  The  temperature  is  developed  after  the  30- 
day  run  (Figure  8),  but  still  close  to  the  day  one  temperature  when  no  data 
is  inserted. 

The  effectiveness  of  the  insertion  of  CAT  data  can  be  illustrated  in  Fig- 
ure 9  which  shows  the  global  volume  rms  values  of  the  differences  between 
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this  experiment  without  data  assimilation  and  the  comprehensive  experiment 
with  data  insertion  EXs-  Figure  9  shows  that  at  an  initial  time  (day  1)  the 
differences  are  very  small.  The  differences  are  close  to  zero  at  day  1  because 
the  temperature  fields  of  EXo  and  EX3  are  the  same  at  day  zero.  The  fields 
are  more  different  after  more  days  of  insertion.  The  growth  of  the  global 
rms  value  of  the  differences  must  be  understood  as  a  "good  sign"  of  the  data 
insertion  process.  Here,  the  inserted  data  acts  as  an  initial  convergence  fac- 
tor to  force  the  model  to  "lose"  the  initial  state  quickly  and  converge  to  the 
"true  ocean"  (see  the  experiments  EXi ,  EXs  and  EX3).  The  global  rms  dif- 
ferences grow  very  fast  in  the  first  ten  days,  slow  down  in  the  next  ten  days, 
and  becomes  stable  at  about  day  30.  A  measure  of  the  time  scale  for  loss  of 
predictability  of  the  model  is  Timepredictabiuty  =~  30  days  in  the  experiment 
EX3  and  other  experiments  (not  shown). 

4.3     The  NE-SW  Assimilation  Experiment  EXi 

The  NE-SW  CAT  slice  section  can  be  seen  in  Figure  10  by  a  chain  of  11 
squares  which  has  a  distance  of  300-350  km.  The  model  is  spun  up  for  30 
days  with  temperature,  salinity,  and  wind  climatologies  which  are  described 
above.  Then  the  EXi  experiment  is  performed  by  the  insertion  of  CAT 
data  for  30  days  in  the  NE-SW  slice  at  one  day  time  step  and  at  all  levels. 
Figures  10-12  and  13-15  show  the  EXi  temperature  fields  at  the  surface  (level 
1),  mixed  layer  level  (level  6),  and  the  bottom  level  (level  15)  for  days  one 
and  30,  respectively.  Comparing  Figure  1  of  the  "true  ocean"  and  Figure  10 
of  the  NE-SW  slice  assimilated  at  the  surface  after  one  day  of  insertion,  it 
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is  seen  that  the  two  surface  temperature  fields  are  very  different  from  each 
other,  both  in  shallow  and  deep  water.  Here,  the  model  accepts  the  CAT 
NE-SW  data  slice  and  spreads  it  out  in  both  directions  but  more  towards  the 
east.  Figure  13  shows  the  field  after  30  days  of  assimilation.  Figure  13  also 
shows  the  spreading  is  increased  in  the  shallow  water  which  may  result  from 
stronger  mixing  and  horizontal  shear  and  its  interaction  with  topography. 

— # 

The  stronger  horizontal  shear  increases  the  horizontal  diffusion  [terms  F  and 
Fe,  in  equations  (21)  and  (24)].  This  is  made  clearer  by  comparing  Figure 
14  of  day  30  assimilation  at  level  6  with  Figure  2  of  the  "true  ocean"  at 
the  same  level.  Here,  the  CAT  data  from  the  NE-SW  assimilated  slice  is 
accepted  by  the  model  and  is  spread  out  towards  the  southeast  where  the 
Gulf  has  shallower  average  depths.  In  the  bottom  layer  the  Figure  6  of  the 
"true  ocean"  and  the  Figure  15  of  the  NE-SW  assimilation  are  very  much 
the  same. 

The  convergence  of  the  NE-SW  CAT  slice  assimilation  experiment  to- 
wards the  "true  ocean"  can  be  measured  by  the  global  rms-error  over  a 
volume  of  the  entire  domain  studied.  Figure  16  shows  the  rms-error  between 
the  "true  ocean"  and  the  NE-SW  slice  experiment  EXi  for  30  days  of  as- 
similation. The  rms-errors  are  0.95,  0.925,  and  0.64  for  days  0,  1,  and  30, 
respectively.  This  means  that  the  rms-errors  of  the  NE-SW  slice  assimilation 
decrease  3  and  33  %  after  one  day  and  30  days  of  assimilation.  The  model 
and  assimilation  scheme  work  very  well  in  this  case.  The  rms-error  curve  is 
almost  linear  and  it  is  expected  that  the  convergence  to  the  "true  ocean" 
would  be  even  better  if  the  insertion  were  made  more  often  and  at  higher 
horizontal  and  vertical  spatial  resolution. 
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4.4      The  NW-SE  Assimilation  Experiment  EXz 

This  experiment  is  very  similar  to  the  NE-SW  CAT  shce  above.  The  NW-SE 
CAT  slice  can  be  seen  in  Figure  17  by  a  chain  of  11  squares  over  a  distance 
of  300-350  km.  In  this  experiment,  the  model  is  spun  up  for  30  days  with 
the  Levitus  temperature  and  salinity,  and  the  ECMWF  wind  climatologies. 
Then  the  model  is  run  for  30  days  with  the  CAT  data  insertion  along  the 
NW-SE  slice  at  one  day  time  interval  for  all  levels  (16  level)  to  perform  the 
assimilation  experiment  EX2. 

Figures  17-19  and  20-22  show  the  EXg  temperature  fields  at  the  level  1 
(surface),  level  6  (mixed  layer),  and  level  15  (bottom  layer)  for  days  1  and 
30,  respectively.  In  general,  it  is  seen  from  these  figures  that  temperature 
fields  are  developing  towards  the  "true  ocean".  As  in  the  experiment  EXi, 
comparing  Figure  1  of  the  "true  ocean"  and  Figure  17  of  the  NW-SE  slice 
assimilated  at  the  surface,  shows  that  the  two  surface  temperature  fields 
are  very  different  from  each  other  everywhere  in  the  region.  The  surface 
temperature  field  after  30  days  of  CAT  assimilation  is  shown  in  Figure  20.  It 
is  seen  that  the  model  accepts  the  CAT  slice  data  and  spreads  it  out  along 
both  sides  of  the  NW-SE  section  at  the  surface.  It  is  clear  from  Figure  5 
of  the  "true  ocean".  Figure  18  of  day  one,  and  Figure  21  of  day  30  that 
at  level  6  the  temperature  of  the  NW-SE  slice  is  shifted  and  more  spread 
out  northwards  in  the  mixed  layer.  As  mentioned  in  the  EXj ,  this  may  be 
the  result  of  the  strong  horizontal  shear  in  the  NW-SE  slice  region  which 
increases  the  horizontal  diffusion  [terms  F  and  F^,  in  equations  (21)  and 
(24)].  In  the  bottom  layer  the  temperature  fields  of  the  "true  ocean"  at  day 
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30  in  Figure  6  and  of  the  NW-SE  slice  assimilation  experiment  EXs  at  day 
30  in  Figure  22  are  very  much  the  same. 

The  convergence  of  the  temperature  field  to  the  "true  ocean"  in  the  NW- 
SW  CAT  slice  experiment  EX2  can  be  measured  by  rms-errors  in  Figure  23. 
This  figure  shows  the  global  rms-error  over  the  volume  of  the  entire  domain 
between  the  "true  ocean"  and  the  the  NW-SE  sHce  experiment  EX 2.  It  is  in- 
teresting to  note  that  the  global  rms-errors  of  EX}  and  EX2  have  nearly  the 
same  distribution  for  30  days  of  assimilation  which  are  0.95,  0.925,  and  0.065 
for  days  0,  1,  and  30,  respectively.  As  in  EXj,  the  rms-errors  of  the  NW-SE 
experiment  decrease  3  and  about  33  %  after  one  day  and  30  days  of  assimila- 
tion, respectively.  One  of  the  reasons  for  the  similarity  in  the  two  rms-error 
distributions  is  that  two  slices  (NE-SW  and  NW-SE)  have  the  same  length 
(300-350  km)  and  are  very  symmetric.  As  in  the  experiment  EXj ,  the  model 
and  the  assimilation  scheme  work  very  well  together  in  this  case.  The  model 
is  greatly  improved  (33  %)  after  30  days  of  CAT  insertion.  It  is  also  seen  that 
the  behavior  of  the  global  rms-errors  is  nearly  insensitive  to  the  change  in  the 
CAT  data  section  location  in  energetic  regions  of  the  Gulf  in  the  experiments 
EXi  and  EX2-  These  two  slices  have  almost  the  same  effectiveness  in  improv- 
ing the  model  errors  on  the  global  space  scale.  These  experiments  also  show 
that  the  internal  forcing  provided  by  the  single  slice  of  inserted  CAT  data  is 
not  localized  in  physical  space  in  inducing  immediate  correlation  between  the 
fully  nonlinear  primitive  equation  model  and  the  "true  ocean".  This  result 
is  very  different  from  the  results  of  Malanotte-Rizzoli  and  Holland  (MRH, 
1988)  who  found  that  the  single  section  of  inserted  data  is  too  localized  in 
physical  space  to  induce  any  immediate  correlation  between  a  QG  model  and 
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the  reference  ocean.  The  reason  for  these  different  results  is  the  differences 
not  only  in  data  and  models  used,  but  also  in  our  smaller  coastal  domain. 
The  rms-error  curve  is  almost  linear  in  this  case  and  it  is  expected  that  the 
convergence  to  the  "true  ocean"  will  be  better  if  the  temporal  and  spatial 
resolutions  of  the  CAT  insertion  is  higher. 

4.5     The  Combined  Slices  Assimilation  Experiment  EX3 

The  more  comprehensive  experiment  is  EXs  which  is  more  realistic  for  coastal 
acoustic  tomography.  The  combined  slices  in  experiment  EX3  include  the 
EXi  slice  and  EXs  slice  together  and  three  other  CAT  slices.  All  these  slices 
together  make  a  parallelogram  without  the  southern  side.  The  EXi  and  EXs 
slices  are  two  diagonals  of  this  CAT  parallelogram  which  is  seen  in  Figure  24. 
The  lengths  of  the  western  and  eastern  sides  of  the  parallelogram  (which  will 
be  called  LON  LEFT  and  LON  RIGHT)  are  about  200  km,  and  the  length 
of  the  northern  side  (which  will  be  called  LAT  slice)  is  about  250  km.  As  in 
two  previous  experiments,  the  model  is  spun  up  for  30  days  with  the  Levitus 
temperature,  salinity,  and  ECMWF  wind  climatologies.  Then  the  CAT  data 
is  inserted  into  the  model,  as  described  in  the  scheme  (31)  and  (32),  for  30 
days  at  one  day  time  intervals  and  at  16  levels  to  perform  the  combined  slices 
assimilation  experiment  EX3. 

The  temperature  fields  of  the  combined  CAT  slices  after  one  day  of  as- 
similation at  the  surface  (level  1),  the  mixed  layer  level  (level  6),  and  the 
bottom  levels  (level  15)  are  shown  in  Figures  24-26.  Figures  27-29  show  the 
temperature  fields  of  the  above  levels  for  day  30.  It  is  seen  from  these  figures 
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that  the  temperature  fields  of  all  levels  at  day  one  are  very  "far"  from  the 
"true  ocean".  As  in  two  other  experiment,  the  temperature  fields  in  this 
experiment  of  levels  1  and  5  are  closer  to  each  other  than  to  the  fields  of  the 
bottom  layer.  This  is  a  result  of  the  surface  layer  and  the  mixed  layer  levels 
being  much  closer  with  stronger  mixing,  especially  in  shallow  water.  After 
30  days  of  the  combined  slices  assimilation,  the  temperature  at  the  surface 
(Figure  27)  is  spread  out  in  all  directions,  but  stronger  in  the  southeastern 
corner  and  in  the  southern  part  of  the  region  where  the  water  is  shallower. 
The  mixing,  horizontal  diffusion,  nonlinearity,  and  interaction  with  topogra- 
phy play  an  important  role  in  the  advection  process  of  temperatures.  At  day 
30  the  SST  field  of  this  experiment  (Figure  27)  converges  well  to  the  "true 
ocean"  of  Figure  4.  It  is  seen  in  all  three  experiments  that  the  temperatures 
of  the  mixed  layer  level  (Figures  14,  21  and  28)  at  day  30  converge  very  well 
to  the  "true  ocean"  in  comparison  with  the  surface  levels.  This  results  from 
the  fact  that  the  mixing  is  strongest  inside  the  mixed  layer  where  the  vertical 
turbulence  has  the  largest  activity.  The  turbulence  exchange  coefficients  Krn 
and  Kh.  in  (21)  and  (24)  have  maximum  values  inside  mixed  layers  and  min- 
imum values  at  the  surface  (see  Ly,  1991).  Inside  the  mixed  layer  both  "top 
down"  and  "bottom  up"  processes  must  be  active  while  only  the  "bottom 
up"  process  is  active  at  the  surface.  These  processes  are  important  in  our 
cases  because  of  the  CAT  data  insertion  by  vertical  slices.  The  temperature 
field  at  level  15  after  30  days  of  assimilation  in  Figure  29  is  almost  identical 
with  the  "true  ocean"  in  Figure  6. 

The  convergence  temperatures  can  be  seen  in  Figure  30  which  shows  the 
global  rms-error  over  a  volume  of  the  entire  domain  between  the  "true  ocean" 
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and  the  the  temperature  in  this  experiment  EXs-  As  expected,  the  conver- 
gence of  the  temperature  field  to  the  "true  ocean"  is  best  in  this  experiment 
in  comparison  with  the  two  experiments  above.  The  rms-errors  drop  dra- 
matically in  the  first  ten  days  of  assimilation  from  0.92  to  0.72  deg  C.  This 
makes  the  convergence  to  the  "true  ocean"  of  22  %  after  the  first  ten  days  of 
assimilation.  The  rms-errors  drop  slower  in  the  next  10  days  from  0.72  at  day 
10  to  0.59  at  day  20.  The  convergence  in  this  second  ten  days  of  assimilation 
is  18  %.  In  the  last  ten  days  of  assimilation  the  rms-errors  drop  slowest  from 
0.59  to  0.50  at  day  30  which  makes  the  convergence  of  15  %.  Totally,  the 
global  rms-errors  drop  is  about  50  %  after  30  days  of  data  insertion.  The 
model  and  the  assimilation  scheme  work  very  well  together.  This  shows  that 
inserting  the  "right"  data  into  a  "wrong"  model  tends  to  drive  the  "wrong" 
model  closer  to  the  "truth"  and  to  reduce  global  rms-errors.  The  inserted 
CAT  data  acts  as  an  initial  convergence  maker  in  the  model.  This  conver- 
gence grows  very  fast.  After  30  days  of  insertion  the  convergence  is  increased 
by  almost  50  %.  This  also  shows  that  the  "true  ocean"  and  the  model  with 
CAT  data  insertion  are  strongly  correlated.  As  in  the  two  above  experiments, 
the  rms-error  curve  is  very  close  to  linear  and  it  is  expected  the  convergence 
to  the  "true  ocean"  will  be  better  if  the  temporal  and  spatial  resolution  of 
the  CAT  insertion  is  higher.  It  is  interesting  to  note  that  convergence  are 
not  proportional  to  the  number  of  CAT  slice  sections  in  an  experiment.  The 
experiment  EXj  has  a  convergence  of  33  %;  EXs  has  a  convergence  about 
33  %,  but  EXs  has  about  50  %). 

It  is  interesting  to  see  the  convergence  of  this  experiment  EXs  to  the 
"true  ocean"  after  another  30  days  of  data  insertion.    Figure  31  shows  the 
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temperature  field  at  the  mixed  layer  level  (level  6)  at  day  60  in  comparison 
with  the  "true  ocean"  at  this  day.  The  temperature  converges  significantly 
to  the  "true  ocean"  at  day  60.  The  temperature  fields  at  this  level  are  almost 
identical.  The  rms-error  distribution  for  days  between  31  and  60  is  shown 
in  Figure  32.  The  rms-errors  between  the  "true  ocean"  and  the  temperature 
field  of  this  experiment  are  reduced  by  17  %  after  30  days  (31-60)  of  insertion. 
Totally,  after  60  days  of  CAT  data  insertion  the  rms-errors  are  reduced  by 
66  %.  The  model  is  improved  by  66  %  after  60  days  of  assimilation. 

Three  other  assimilation  experiments  with  CAT  insertions  along  the  north- 
ern (LAT),  western  (LON  LEFT)  and  eastern  (LON  RIGHT)  sides  of  the 
above  parallelogram  (Figure  24)  have  been  carried  out  to  study  the  effective- 
ness of  various  local  slices.  The  model  spreads  these  CAT  slices  out  differ- 
ently, and  this  strongly  depends  on  the  horizontal  diffusion,  mixing  process, 
and  nonlinear  interaction  with  topography.  At  the  surface,  after  30  days  of 
assimilation  the  LAT  slice  is  spread  out  along  both  sides,  but  more  towards 
the  north  (not  shown).  The  LON  LEFT  slice  is  spread  out  greater  in  the 
northern  part  towards  the  east,  after  30  days  (not  shown).  And  at  the  sur- 
face the  LON  RIGHT  slice,  after  30  days,  is  spread  out  greater  in  the  middle 
towards  the  west  (not  shown).  Although  these  CAT  slices  have  various  tem- 
perature advection,  but  the  their  global  rms-errors  distributions  are  almost 
the  same  (not  shown).  Their  rms-errors  are  reduced  by  27  %  after  30  days 
of  assimilation.  These  experiments  confirm  that  there  is  no  most  effective 
slice  between  different  CAT  slice  sections.  All  slices  have  almost  the  same 
"weight"  in  forcing  the  coast  model  in  converging  to  the  "true  ocean". 
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5     SUMMARY 

Coastal  acoustic  tomographic  (CAT)  data  is  used  in  "twin"  experiments  by 
direct  data  insertion.  A  series  of  assimilation  experiments  has  been  carried 
out  and  discussed  in  which  data  are  inserted  along  various  CAT  slices  into  a 
coastal  ocean  primitive  equation  model  for  the  Gulf  of  Sirte.  The  model  is 
developed  for  the  whole  Mediterrainean  (MED)  Sea  and  is  used  for  studying 
a  direct  insertion  technique.  The  model  has  realistic  coastlines  and  topogra- 
phy. This  model  also  has  a  coastal-following,  curvilinear,  nearly  orthogonal 
horizontal  coordinate  system,  fully  active  thermodynamics,  a  second  order 
turbulence  closure,  16  bottom-following,  sigma  vertical  levels,  and  an  aver- 
age resolutions  of  about  10  by  10  km.  The  model  was  spun  up  for  30  days 
with  Levitus  temperature  and  salinity  and  ECMWF  wind  climatologies  and 
then  run  for  a  year  to  generate  the  "true  ocean".  In  the  assimilation  experi- 
ments, after  a  30-day  spinup,  the  model  was  run  for  30  days  with  an  insertion 
time  interval  of  one  day  at  all  16  vertical  levels  (except  the  comprehensive 
experiment  EXs  which  was  run  for  60  days). 

The  results  indicate  that  the  model  works  well  with  the  assimilation 
scheme  to  accept  data  slices  and  spread  it  out  nearly  in  all  directions  in 
the  Gulf  of  Sirte.  The  inserted  data  acts  as  an  initial  convergence  maker 
which  forces  the  model  to  "lose"  the  initial  state  quickly  and  converge  to 
the  "true  ocean".  The  global  rms-errors  in  the  studied  domain  between  the 
"true  ocean"  and  the  temperature  fields  of  the  comprehensive  experiment 
EXs  decrease  by  about  50  %  after  30  days  and  66  %  after  60  days  of  as- 
similation. The  rms-errors  of  two  other  experiments  with  diagonal  slices  are 
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reduced  by  33  %,  and  of  all  other  sides  of  the  parallelogram  are  reduced  by 
27  %  after  30  days  of  insertion.  This  shows  that  inserting  the  "right"  data 
into  a  "wrong"  model  tends  to  drive  the  "wrong"  model  closer  to  the  "truth" 
and  to  reduce  global  rms-errors.  The  good  convergence  also  shows  that  the 
"true  ocean"  and  the  model  with  CAT  data  insertion  are  strongly  correlated. 
The  rms-error  curve  is  very  close  to  linear  in  the  30  days  of  insertion.  Then, 
it  is  expected  that  the  convergence  to  the  "true  ocean"  will  be  better  if  the 
temporal  and  spatial  resolutions  of  the  CAT  insertion  are  higher. 

The  temperature  is  spread  out  greater  in  the  shallow  water  which  may 
be  the  result  of  stronger  mixing  and  horizontal  shear  and  its  interaction  with 
topography,  which  increases  horizontal  diffusion.  It  is  seen  in  all  experiments 
that  the  temperatures  of  the  mixed  layer  level  at  day  30  converge  very  well 
to  the  "true  ocean"  in  comparison  with  the  surface  levels.  This  results  from 
the  fact  that  the  mixing  is  strongest  inside  mixed  layers.  The  turbulence 
activities  are  strongest  inside  mixed  layers  and  weakest  at  the  surface.  Inside 
mixed  layers  both  "top  down"  and  "bottom  up"  processes  must  be  active 
while  only  the  "bottom  up"  process  is  active  at  the  surface.  These  processes 
are  important  in  CAT  data  assimilation  because  of  insertions  by  vertical 
slices. 

The  behavior  of  the  global  volume  rms-errors  is  nearly  insensitive  to  the 
change  of  the  CAT  data  slice  location  in  the  more  energetic  regions  of  the 
Gulf.  Although  side-slices  and  diagonal-slices  of  the  parallelogram  have  var- 
ious temperature  advection,  their  global  rms-errors  distributions  are  almost 
the  same.  Their  rms-errors  are  reduced  by  27  %  and  33  %  after  30  days  of 
assimilation.  There  is  no  the  most  effective  slice  between  different  local  CAT 
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slices.  All  slices  have  almost  the  same  "weight"  (effectiveness)  in  improving 
the  model  errors  on  the  global  space  scale  of  the  Gulf  of  Sirte.  The  conver- 
gence is  not  proportional  to  the  number  of  CAT  slices  in  an  experiment. 

The  experiments  also  show  that  the  internal  forcing  provided  by  the  single 
slice  of  inserted  CAT  data  induces  an  immediate  correlation  between  the  fully 
nonlinear  primitive  equation  model  and  the  "true  ocean"  that  is  not  localized 
in  physical  space.  This  result  is  very  different  from  that  of  Malanotte-Rizzoli 
and  Holland  who  found  that  the  single  section  of  inserted  data  is  too  local- 
ized in  physical  space  to  induce  any  immediate  correlation  between  a  QG 
model  and  the  reference  ocean  (MRH,  1988).  The  reason  for  our  different 
conclusions  is  due  to  the  data  and  models  used  and  to  our  smaller  coastal 
domain. 

Overall,  a  direct  insertion  scheme  for  CAT  data  works  very  well  with  a 
coastal  multilevel,  primitive  equation  model  with  complete  thermodynamics, 
a  second  order  turbulence  closure,  and  realistic  coastlines  and  topography. 
A  comparison  of  this  scheme  with  a  modified  Kalman  filter  using  CAT  is 
recomended  for  future  investigation. 
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7     Figure  Captions 

Table  1  Block  diagram  of  the  simulation  twin  experiments  for  studying 
cocistal  acoustic  tomography  data  assimilation. 

Fig.  1  Temperature  field  of  the  "true  ocean"  at  the  sea  surface  (level  1)  on 
day  1. 

Fig.  2  Same  as  figure  1,  but  for  level  6  (a  mixed  layer  level). 

Fig.  3  Temperature  field  of  the  "true  ocean"  on  day  1  for  level  15  (the 
bottom  layer). 

Fig.  4  Temperature  of  the  "true  ocean"  at  the  sea  surface  (level  1)  on  day 
30. 

Fig.  5  Same  as  Fig.  4  but  for  level  6. 

Fig.  6  The  "true  ocean"  at  level  15  (the  bottom  layer)  on  day  30. 

Fig.  7  Temperature  of  the  sea  surface  for  day  1  in  the  experiment  EXg  with 
no  data  assimilated. 

Fig.  8  Same  as  Fig.  7  but  for  day  30. 

Fig.  9  Root-mean-square  (rms)  values  of  the  differences  between  the  exper- 
iment with  the  insertion  of  five  CAT  slices  {EXs)  and  the  experiment 
without  data  insertion  {EXq))- 
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Fig.  10  Temperature  at  the  sea  surface  on  day  1  of  the  experiment  EXj 
with  a  single  NE-SW  tomographic  section.  The  tomographic  section 
can  be  seen  as  a  chain  of  11  squares. 

Fig.   11   Same  as  Fig.  10  but  for  level  6  (a  mixed  layer  level). 

Fig.   12  Same  as  Fig.  10  but  for  level  15  (bottom  layer) 

Fig.  13  Temperature  at  the  sea  surface  for  day  30  of  the  experiment  EXi 
with  a  single  NE-SW  tomographic  section. 

Fig.  14  Same  as  Fig.  13  but  for  level  6. 

Fig.   15  Same  as  Fig.  13  but  for  level  15. 

Fig.  16  Global  volume  rms  error  versus  time  for  a  period  of  30  days  be- 
tween the  nowcast  temperature  fields  in  EXi  and  the  "true  ocean" 
temperature  fields. 

Fig.  17  Temperature  at  the  sea  surface  (level  1)  on  day  1  of  the  experiment 
EX2  with  a  single  NW-SE  tomographic  section.  The  tomographic  sec- 
tion can  be  seen  as  a  chain  of  11  squares. 

Fig.  18  Same  as  Fig.  17  but  for  level  6  (a  mixed  layer  level). 

Fig.   19  Same  as  Fig.  17  but  for  level  15  (bottom  layer). 

Fig.  20  Temperature  at  the  sea  surface  on  day  30  in  EX2  with  a  NW-SE 
tomographic  section. 

Fig.  21  Same  as  Fig.  20  but  for  level  6. 
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Fig.  22  Same  as  Fig.  20  but  for  level  15. 

Fig.  23  Global  volume  rms  error  versus  time  for  a  period  of  30  days  be- 
tween the  nowcast  temperature  fields  in  EX 2  and  the  "true  ocean" 
temperature  fields. 

Fig.  24  Day-1  sea  surface  temperature  of  the  comprehensive  experiment 
EXs  with  five  tomographic  sections  forming  a  parallelogram  enclosing 
a  coastal  region.  The  parallelogram  (without  a  southern  section)  can 
be  seen  as  five  chains  of  brightened  squares. 

Fig.  25  Same  as  Fig.  24  but  for  level  6. 

Fig.  26  Same  as  Fig.  24  but  for  level  15. 

Fig.  27  Sea  surface  temperature  on  day  30  in  the  experiment  EX3  with  the 
assimilation  of  five  tomographic  sections. 

Fig.  28  Same  as  Fig.  27  but  for  level  6. 

Fig.  29  Same  as  Fig.  27  but  for  level  15. 

Fig.  30  Global  volume  rms  error  versus  time  for  a  period  of  30  days  be- 
tween the  nowcast  temperature  fields  in  EX3  and  the  "true  ocean" 
temperature  fields. 

Fig.  31  Temperature  at  a  mixed  layer  level  (level  6)  in  the  comprehensive 
experiment  EXs  right  after  60  days  of  data  insertion.  Temperature 
of  the  "true  ocean"  for  the  same  time  (i.e.,  day  60)  and  same  level  is 
shown  in  the  lower  panel. 
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Fig.  32  Global  volume  rms  error  versus  time  from  day  31  to  day  60  between 
the  nowcast  temperature  fields  in  EXs  and  the  "true  ocean"  tempera- 
ture fields. 
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Day  =  60  Level  =  6  True  Ocean 
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Fig.   31 


Vertical  tomographic  combined  slices  assimilated 
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